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1. The determination of the AS = 1, CP-violating parameter e'/e in neutral 
kaon decays has attracted considerable experimental and theoretical interest (for a 
complete list of references see, for instance, ref. The possibility of drastic can- 
cellations among QCD- and electroweak-induced operators ||, occurring for large 
values of the top quark mass, has spurred new and more accurate theoretical investi- 
gations [[| |]. On the experimental side, the present sensitivity of 1 part in 10 3 is 
not conclusive and an effort in improving it by one order of magnitude is under way. 

The theoretical framework for the study of e'/e is the effective field theory. The 
local Hamiltonian for AS" = 1 transitions can be written, for /i < m c , as |1 

G 

n = ^KJ2[ z i(^) + TVi(t4\Qi(p) ■ (!) 

The list of the effective operators Qi {i = \ — 10) is reported in refs. [[| |J], whose 
notation we follow closely and where the reader may find a complete discussion of 
the basic tools we use in the present analysis. We just recall that Q 1>2 stand for 
the "tree-level" PF-induced current-current operators, Q3-6 for the QCD penguin 
operators and Q7-10 for the electroweak penguin (and box) ones. 

The Wilson coefficients 24,2 (/■*) run from mw to m c via the corresponding 2x2 
sub-block of the 10 x 10 anomalous dimension matrix, while Zj(/i) = for i = 3 — 10. 
From fi = m c down, as the charm-induced penguins come into play, ^(/i) evolve, 
given the proper matching conditions, with the full anomalous dimension matrix. 
The Wilson coefficients Vi(fi) (y«(/i) = fj(/i) — Zi(fJ>)) arise at mw due to integration 
of the W and top quark fields. They coincide with Zi(fi) for i — 1,2, the information 
about the top quark being encoded in the i = 3 — 10 components. Finally, Aj = V ic [V* s , 
where V is the Kobayashi-Maskawa (KM) matrix, and r = —X t /X u . 

Let us mention that, according to standard conventions, Im A„ = and the short 
distance component of e'/e is thus determined by the Wilson coefficients j/j. Notice 
that, following the approach of ref. 0, yi(fi) = ?/2(a0 = and the effect of C/1,2 
appears only through the linearly dependent operators Q/i^io- 

In this letter, we would like to discuss the possible relevance of two extra operators 
that have so far been neglected, namely 



<2u = Sh m * s^ut a G^(l - 75) d (2) 



and 



eQ 



Q12 = ^rn s sa^F^il -j 5 )d, (3) 
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which account for the magnetic and electric dipole part of, respectively, the QCD 
and electromagnetic penguin operators. In eq. (||), t a are the SU(3) C generators 
normalized as Tr(t a t b ) = {l/2)5 ab , while in eq. (|3|) Qd = —1/3 is the charge of the 
down quarks; = (V 2 )bv, lA- 

These operators have been left out in the past, together with Z penguins and 
electroweak box diagrams, because their Wilson coefficients exhibit, for light quark 
masses, a power-like GIM suppression 0, to be compared with the leading soft loga- 
rithmic behavior of the "monopole" component of the gluonic and photonic penguins. 
However, for a heavy top quark, there is no reason a priori to discard these contribu- 
tions, which should be considered on the same grounds as Z penguins and electroweak 
box diagrams. In addition, Qu and Q12 receive a large QCD renormalization due to 
the mixings with the other operators, as has been shown in the case of b — > 57 |7|] . 

The potential relevance of the dipole components of the gluonic and photonic 
penguins for e' / e has been emphasized by one of us a few years ago || . Only recently 
has a quantitative attempt to estimate the relevance of Qu appeared ||. However, 
in ref. |§ , only the multiplicative renormalization — which is by no means the leading 
effect — is taken into account. We also disagree on the evaluation of the hadronic 
matrix element that, we believe, is overestimated in ref. || by more than one order 
of magnitude (by a factor 8 using their method). 

The evaluation of the matrix elements is in fact the crucial point in determining 
the relevance of the two additional operators. As we shall see, while the contribution 
of the electromagnetic penguin Q 12 seems to be vanishingly small, that of the gluonic 
Qu should not be neglected. In fact, its Wilson coefficient turns out to be sizeably 
enhanced by mixing and the matrix element is not dramatically suppressed with 
respect to the other ten operators. 

2. In order to compare meaningfully the relevance of the new contributions with 
the traditional ones, from now on we strictly follow the analysis presented in ref. 0|. 
The parameter e'/e can be obtained by means of the Hamiltonian in eq. (jl]), which 
yields § 
F > 

- = 10- 4 
e 

where 

P (i/2) = r J^Vi (2tt, I ^0\Qi\K°)(l-n v+v/ ) (5) 
pim = Lj2 yi (2ir, I = 2\Q % \K») . (6) 
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Im A t 
1.7 x 10- 4 



p(l/2) _ p(3/2) 



(4) 



We take, as input values for the relevant quantities, the central values given in ap- 
pendix C of ref. ||. This will allow us to reproduce in the ten-operator case the 
central values of the leading order (LO) results, as given in appendix B of ref. ||. In 
particular, we take 

r = 1.7 ^ ~ 594 GeV -3 , lu = 1/22.2 , Vt^ = 0.25 ; (7) 
2 |e| Re Aq 

Im A t is determined from the experimental value of e as an interpolating function of 
m t . For instance, taking the KM phase 6 km in the first quadrant, we find, for the 
central value, 

ImAi = 2.77 x 1(T 4 a;" 0365 , (8) 
where x t = mf/rriyy 

The value of the Wilson coefficients yu and ?/ 12 at the hadronic scale of 1 GeV 
can be found by means of the renormalization group equations. Denoting generically 
by C(fi) the vector of Wilson coefficients, its scale dependence is given by 



9 n , \ 9 




iV.^^.J 1 ,"), (9) 



where /3(g) is the QCD beta function and a the electromagnetic coupling (the running 
of a is being neglected). At the leading order we have 

7( 9 2 ,«) = ^7f» + ^7<°\ (10) 

where 7^ governs the QCD and 7^ the electromagnetic running. 

The 10 x 10 mixing matrix of the anomalous dimension for Q1-10 has recently been 
evaluated at the next-to-leading order in refs. 0, M. On the other hand, the matrix 
of the strong anomalous dimensions of Q11 and Q12, and their QCD- induced mixing 
with Qi_6> can be borrowed from the existing calculations for the b — ► 57 decay J7[ 
(we will use the most recent and complete analysis of ref. flOfl). Assembling the known 
results, we thus write the leading order 12 x 12 anomalous dimension matrix *yf' T as 
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-6/N 6 

6 -6/JV 

-2/3JV -22/3JV 6-2n//3JV 
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2/3 4/3 2n//3 -6 
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-6/JV 










32JVn//9 + 50n//9JV 


x 7 


x 8 


*9 


^10 


AN - 8/N 





) 


-12n/ +8n f /9)C f 




Y 8 


Yg 




8C f 


8C f 



where N is the number of colors, n/ is the number of active flavors, and rif = 
nd + (Q u /Qd)n u , with Q u = 2/3, Q d = -1/3, n u (n d ) being the number of up (down) 
active quarks; finally, Cf = (iV 2 — 1)/2N. 

A few comments are in order. The 10 x 10 part of the matrix ( |Tl~D is identical to 
that used in refs. || || for the leading one-loop order result. The two extra columns 
and rows represent the mixing of the first ten operators with the two new ones, which 
takes place first at the two-loop level. We have taken for these entries the 't Hooft- 
Veltman (HV) scheme results || |J. In this way, no finite additional contributions 
to the renormalization of yn and yu arise at the various quark thresholds (for a 
discussion see ref. HT0|). The entries labelled with X 7 _ w , Y 7 _ 10 represent the mixings 
of, respectively, Qu and Q 12 with the electroweak penguins, which have not yet been 
computed. We set them equal to zero, expecting that the leading contribution to 
the running of the Wilson coefficients yn and yu arises from the mixing with the 
current- current operators and the gluonic penguins. 

Since we want to fully reproduce the leading-order results for the 10 x 10 operator 
system, we have included, in the running of the Wilson coefficients, the effect of 
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the electromagnetic renormalization following the procedure described in ref . |3j . To 
the 10 x 10 electromagnetic anomalous dimension matrix, which can be found in 
refs. |3|, |J, we have simply added two rows and two columns of zeros, thus neglecting 
the electromagnetic running of Qn and Qi 2 - 

The fact that the mixing of Qu and Qu with the other operators arises at the 
two-loop level introduces an explicit scheme dependence in our analysis. From the 
next-to-leading results of refs. |3|, f|, it appears however that the inclusion of next-to- 
leading effects reduces the values of e'/e compared to the leading-order result. Since 
the presence of the two extra operators does not affect the running of 2/i_i (/i) (see 



eqs. (P)- (|TT1) ), we conclude that our simplified leading-order procedure will at most 
underestimate the effective weight in e'/e of the new contributions. 

For the purpose of comparison, we take as initial conditions for fi_io(^iy) the 
leading-order values of ref. 0, which coincide with the HV results when neglecting 
0(a s ) corrections to v\^{mw) (strictly speaking at the leading order one should set 
v 3 _io(mw) = 0). For what concerns the new coefficients v n,i 2 (m\y) we have 

Vn{m w ) = -E'{m 2 t /m 2 w ) and v 12 (m w ) = -D'(m 2 /m 2 w )/Q d , (12) 

where 

Qry»2 ,y. 3 C /y» 2 O ,-y. 

E (x) = — rrlnx : — — (13) 

v ; 2(1 -x) 4 4(1 -x) 3 V ; 

x 2 (2-3x) 8x 3 + 5x 2 - 7x 

D (x) = — ; n-mx ; — — . 14 

v ; 2(1 -x) A 12(1 -xf v ; 

In Table 1 we report the results for 2;i_i 2 (l GeV) and 2/3-12(1 GeV) (remember 
that 2/1,2 (Ai) = 0) compared to their initial values, for m t equal to 130 and 170 GeV and 
J^qcd equal to 200 and 300 MeV. As expected, a comparison between our results and 
those of Table 1 of ref. shows that we fully agree on the values of the renormalized 
coefficients for the first ten operators. Regarding f 11,12 (= 2/11,12 + -211,12), we note 
that the effect of operator mixing induces a renormalization that is a factor of 4-5 
larger than that induced by multiplicative running alone (which roughly reduces by 
a factor of 2 the initial Wilson coefficients). 

In order to ascertain the importance of Qu and Q12 in the estimate of e'/e, we 
now have to address the question of the evaluation of their hadronic matrix elements. 
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A (4) 


0.2 GeV 


0.3 GeV 


m t 


130 GcV 


170 GcV 


130 GcV 


170 GcV 


Z\ 


-0.587 (0.0) 


-0.715 (0.0) 


Z2 


1.319 (1.0) 


1.409 (1.0) 


Z3 


0.004 


0.005 


Zi 


-0.010 


-0.014 


Z5 


0.003 


0.004 


Z6 


-0.011 


-0.015 


z 7 /a 


0.005 


0.008 


z 8 /a 


0.001 


0.001 


zg/a 


0.005 


0.009 


zw/oe 


-0.001 


-0.001 


Zll 


-0.038 


-0.045 


Zl2 


0.431 


0.582 


2/3 


0.027 (-0.0003) 


0.028 (0.0004) 


0.034 (-0.0003) 


0.035 (0.0004) 


2/4 


-0.048 (0.0018) 


-0.049 (0.0014) 


-0.056 (0.0020) 


-0.057 (0.0015) 


2/5 


0.011 (-0.0006) 


0.012 (-0.0005) 


0.012 (-0.0007) 


0.013 (-0.0005) 


2/6 


-0.078 (0.0018) 


-0.080 (0.0014) 


-0.102 (0.0020) 


-0.104 (0.0015) 


ViIol 


-0.025 (0.091) 


0.025 (0.151) 


-0.017 (0.091) 


0.033 (0.151) 


2/s/a 


0.053 (0.0) 


0.109 (0.0) 


0.074 (0.0) 


0.148 (0.0) 


2/9/" 


-1.160 (-0.793) 


-1.554 (-1.094) 


-1.222 (-0.793) 


-1.644 (-1.094) 


yw/ot 


0.488 (0.0) 


0.663 (0.0) 


0.592 (0.0) 


0.806 (0.0) 


2/11 


-0.328 (-0.172) 


-0.340 (-0.193) 


-0.350 (-0.172) 


-0.362 (-0.193) 


2/12 


2.070 (0.965) 


2.161 (1.158) 


2.166 (0.965) 


2.245 (1.158) 



Table 1: Wilson coefficients at [i — 1 GeV, in the leading order, as explained in the 
text (a = 1/128). The corresponding values at /i — m w are given in parenthesis. In 
addition, at fj, — m c we have 2 3 _i 2 (m c ) = 0. 
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3. Most of the uncertainties in the estimate of e'/e arise from the evaluation 
of the hadronic matrix elements. For the operators Qi-io, we follow the strategy 
of ref. where the various matrix elements are evaluated by means of the 1/N 
expansion and soft-meson methods. Overall coefficients B^ 2 ^ and B^ 2 ^ (i = 1 — 10) 
parameterize our level of ignorance of their normalization scale, scale dependence and 
error estimate. Since some of the matrix elements are phenomenologically determined 
using CP-conserving data (AI = 1/2 rule), and further relations among the B^s are 
advocated, the final result is parameterized in terms of two coefficients: B^ 2 ^ and 
B^ 2 \ which take 1 as central value. The inclusion of Qu and Q12 requires three 
additional effective parameters: B^ 2 \ ByJ 2 ^ and B^ 2 ^ ' . 

In order to evaluate (tt + 7t~\ Qu \K°) and (7r + 7r~| Q12 | -K" ) , we use effective chiral 
lagrangian techniques. The lowest order chiral representation for these operators van- 
ishes because of the cancellation between the direct transition and the pole diagram 
induced by a non- vanishing (0| Qu \K°) matrix element fill . The leading next order 
(O (p 4 )) contribution arises from the bosonization of Qu into Lq u , where 

L Qll = G Qll Tr { [^M q X+ + A + EMJ] d^d^} . (15) 

In the formula above A + = (A6 + i^i)/^ is the octet AS = 1 projector, M. q is the 
quark current mass matrix, S = exp (y/2iU/ ' F^j and II = J2a=i * s ^ ne USU& 1 
matrix of chiral SU(3) Goldstone bosons. 

The constant Gq u can be computed by considering a model for the QCD effective 
action at long distances [fL^I ; it is found to be [|13| 

G Qll = --^{0\qq\0) , (16) 

where (0|gg|0) is the quark condensate that we take to be — F]^m 2 K /[2.{rn s + m u )). 
More details on this computation will be given elsewhere [I3|. The final result can 
be written as 

Qn\K«) = --^-^^ml m lB^\ (17) 
lo7r z 2 m s + m u r£ 



7T 7T 



with B[\ /2) = 1. 

The corresponding matrix element for Q12 is very much suppressed because it 
is proportional to the photon condensate and it must therefore be very small, if not 
vanishing; we will neglect it altogether in what follows. Nevertheless, we must bear in 
mind the possibility that the same electromagnetic dipole penguin operator could give 
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a sizeable contribution, for instance, once saturated with an external quark current 
and written as a four-fermion operator. Such a possibility is under investigation. 

Notice in eq. (|D]) the presence of the factor 1/167T 2 , reminding us that the leading- 
order mixing between Qi-io an d Qii appears at the two-loop level. 

Since (tt+tt- \Qn\K } = (tt 7v \Q u \K } we finally obtain 

(2ti, I = 0\Q u \K°) = ^(vr+vr-IQnlK ) 

(2tt, I = 2\Q u \K°) = 0, (18) 

which enter in the determination of e'/e (see eq. (|j)). 

We feel that our estimates for the hadronic matrix elements of the operator Qu 
and Q12 are consistent with those in ref. || for the other ten operators. 

4. We can now discuss the effect of the various operators in determining the size 
of e'/e, and identify the role of the new contributions. We follow here the analysis of 
ref. 0, where the hadronic matrix elements are evaluated using the 1/N expansion 
and soft-meson methods. As previously mentioned, the final results are written in 
terms of effective coefficients B^' 2 ^ and B^' 2 ' , which encompass our lack of knowledge 
on scale normalization, scale dependence and goodness of the method. The matrix 
elements of Q\ and Q2 can however be determined phenomenologically from the 
experimental values of Re/Lj and ReA 2 , so as to reproduce the AI = 1/2 rule. In 
particular, in ref. |1 it is found that -B^lo ~ 5-8, which is about three times larger 
that the 1/N result. Related to this coefficient is the value of B^f^,, which we find 

equal to 19.3 and 13.5 for Js}q CD = 0.2 and 0.3 GeV, respectively. Correspondingly, 

(3/2) 

B\°£% = 0.48 and 0.50. Relations among the other coefficients are advocated in 
ref. 0, depending on the relevance and the role of the various operators, so as to 
reduce, in the ten-operator case, the description of the hadronic sector to two effective 
parameters: B^ 2 ^ and B^ 2 \ whose 1/N value is 1. 

Since the determination of B\ and B2 is best achieved at fi = m c 0, all the 
hadronic matrix elements are assumed to be evaluated at that scale and renormalized 

down to 1 GeV via their anomalous dimension matrix. 

(1/2) (1/2) (3/2) 

We proceed analogously by setting B\{ = 1 and B\ 2 = B l2 = at jj = m c 

and using the 12 x 12 QCD and electromagnetic evolution matrices to evolve all the 
hadronic matrix elements to the 1 GeV scale. Since the anomalous dimension matri- 
ces, which govern the evolution of the hadronic matrix elements, are the transpose of 
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those evolving the Wilson coefficients, we now find that the presence of (3n,i2 affects 
the renormalization of the first ten operators. On the other hand, the evolution of 
Qn,i2 is determined solely by their 2x2 anomalous dimension matrix, which imply 
that the matrix element of Q12 remains vanishing. 

As a consequence of the previous remarks, our results for the individual contribu- 
tions of the operators Qi-io t° e V e differ (slightly) from those reported in ref. |J. We 
have however checked that, in the ten-operator case, we reproduce their leading-order 
results exactly. 

We have chosen to illustrate numerically our conclusions as tables. Tables 2, 3 
and 4 show the contributions to e'/e of each operator, for different choices of Aqcd, 

(\/2) (3/2) 

Bq and .Bg . The first ten contributions are also partially grouped in a "positive" 
gluonic component versus a "negative" electroweak component, which shows the "su- 
perweak" behavior of e'/ e within the standard model, as the top mass increases. The 
total effect in the twelve-operator case is then compared with the ten-operator case 
result (we have filled in those data that are not explicitly available from ref. H). 

These tables suggests that Qn should be consistently included in any estimate of 
e'/e. Its sign is the same as that of Qq and therefore makes e'/e larger; for large m t , 
the "super-weak" regime is accordingly shifted to higher values. 

Let us conclude by mentioning the importance of completing the leading-order 
(two-loop) calculation of the anomalous dimension matrices — which is important also 
for the AB = 1 processes. Moreover, a consistent discussion of the short-distance 
part of the present analysis should include the next-to-leading order. This however 
implies computing the mixings of Qn } i2 with the other operators at the three-loop 
level — a truly formidable task. 

We thank J.O. Eeg, E. Franco, S. Narison, P. Nason, N. Paver, A. Pich and 
Riazuddin for helpful discussions. 
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e'/e x 10 4 (Leading Order) 


m t 


130 GeV 


170 GeV 


200 GeV 


230 GeV 


Q 3 


0.5 


9.4 


0.4 


7.8 


0.4 


7.2 


0.3 


6.4 


Q 4 


-5.2 


-4.4 


-3.9 


-3.6 


Q 5 


-0.6 


-0.5 


-0.4 


-0.4 


Qe 


14.7 


12.3 


11.0 


10.1 


Qi 


0.2 


-1.5 


-0.4 


-4.1 


-0.7 


-6.0 


-1.1 


-7.9 


Q 8 


-3.6 


-5.9 


-7.6 


-9.2 


Q9 


3.0 


3.4 


3.6 


3.8 


Qio 


-1.1 


-1.2 


-1.3 


-1.4 


Qn 


0.5 


0.4 


0.3 


0.3 


Qi — Q12 


8.4 


4.2 


1.4 


-1.2 


Ref. f| 


8.0 


3.8 


1.0 


-1.5 



Table 2: Anatomy of e'/e for A ( q ] cd = 300 MeV, B$ /2) = B^ /2) = B{\ /2) = 1. 
The contribution of each operator is shown at /i = 1 GeV, together with partial 
grouping of the gluonic and electroweak sectors. The contribution of Qu is being 
neglected. The total effect is compared with the corresponding leading-order result 
of the ten-operator case (last line). 
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e'/e x 10 4 (Leading Order) 


m t 


130 GeV 


170 GeV 


200 GeV 


230 GeV 


Q 3 


0.4 


6.7 


0.3 


5.6 


0.3 


5.1 


0.3 


4.6 


Qa 


-5.1 


-4.2 


-3.8 


-3.5 


Q 5 


-0.6 


-0.5 


-0.4 


-0.4 


Q 6 


12.0 


10.0 


9.0 


8.2 


Qi 


0.4 


-0.3 


-0.3 


-2.7 


-0.7 


-4.2 


-1.2 


-5.9 


Q 8 


-2.7 


-4.6 


-5.9 


-7.3 


Q9 


2.9 


3.2 


3.5 


3.7 


Qio 


-0.9 


-1.0 


-1.1 


-1.1 


Qn 


0.4 


0.4 


0.3 


0.3 


Qi — Q12 


6.8 


3.3 


1.0 


-1.1 


Ref. § 


6.4 


3.0 


0.7 


-1.4 



Table 3: Same as in Table 2 for A%> CD = 200 MeV, Bq ' 1 = B%' ' = 1. 
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e'/e x 10 4 (Leading Order) 


m t 


130 GeV 


170 GeV 


200 GeV 


230 GeV 


Q 3 


0.5 


5.8 


0.4 


4.8 


0.4 


4.1 


0.4 


3.9 


Qa 


-5.2 


-4.4 


-4.0 


-3.6 


Q 5 


-0.4 


-0.4 


-0.3 


-0.3 


Q 6 


11.0 


9.2 


8.2 


7.5 


Qi 


0.2 


-0.5 


-0.3 


-2.4 


-0.6 


-3.6 


-0.8 


-5.1 


Q 8 


-2.6 


-4.3 


-5.6 


-6.8 


Q9 


3.0 


3.4 


3.6 


3.9 


Qio 


-1.1 


-1.2 


-1.3 


-1.4 


Qn 


0.5 


0.4 


0.3 


0.3 


Qi — Q12 


5.7 


2.8 


0.8 


-1.0 


Ref. § 


5.3 


2.4 


0.5 


-1.3 



Table 4: Same as in Table 2 for A ( q> cd = 300 MeV, Eg'** = 5^ =0.75. 
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